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February  14,  1995 

Consequences  of  Kinetic  Dispersion  on  the 
Electrochemistry  of  an  Adsorbed  Redox-Active  Monolayer 


Gary  K.  Rowe,  Michael  T.  Carter,  John  N.  Richardson,  Royce  W.  Murray* 

Kenan  Laboratories  of  Chemistry,  Univ.  of  North  Carolina,  Chapel  Hill,  NC  27599-3290 


ABSTRACT 

Mixed  monolayers  of  (ferrocenylcarboxy)-alkanethiol/n-alkanethiol  have  been 
investigated  electrochemically  in  2:1  (v:v)  chloroethane:butyronitrile  solvent  in  the 
temperature  range  of  120K  to  150K.  Cyclic  voltammetry  of  these  monolayers  shows  large 
oxidation-reduction  peak  potential  separations  indicative  of  electron  transfer  rate  control. 

The  voltammetric  waveshapes  are  also  broadened;  this  and  curved  log[i]  vs.  time  transients 
observed  in  potential  step  experiments  are  interpreted  as  a  dispersion  in  the  reaction  rates  of 
the  ferrocene  sites.  This  paper  considers  origins  and  three  models  for  such  kinetic 
dispersion:  (i)  Using  simulations,  the  observed  kinetic  dispersion  effects  can  be  successfully 
represented  by  a  Gaussian  distribution  among  the  formal  potentials  E°’  of  the  surface  redox 
sites.  While  only  an  apparent  kinetic  dispersion  (having  a  thermodynamic  origin),  we  show 
by  simulations  that  its  presence  affects  potential  step  log[kApp.,]  vs.  t]  plots,  depressing  the 
apparent  reorganizational  barrier  energies  (X)  and  elevating  the  apparent  rate  constants  (k°), 
consistent  with  previous  experimental  observations.  Similarly,  cyclic  voltammetric 
simulations  with  a  Gaussian  E°’  distribution  give  excellent  fits  to  experimental 


voltammograms  with  mid-point  average  rates  (that  with  voltammograms  can  be  simulated  to 
fit  both  the  experimental  waveshape  and  AEpeak)  that  are  roughly  6-fold  smaller  than  the 
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average  rate  (determined  from  a  fit  to  the  experimental  AEpeak  assuming  a  homogeneous 
population).  The  temperature  and  chain  length  dependence  of  simulations  are  also  consistent 
with  experimental  observations  and  indicate  that  the  dispersion  has  little  effect  on  accurate 
determination  of  X  (from  an  activation  analysis)  or  ^  (from  a  plot  of  log(k°)  vs.  chain  length), 
(ii)  A  Gaussian  distribution  of  reorganizational  energies,  which  is  a  real  kinetic  dispersion, 
has  consequences  on  the  appearance  and  the  analysis  of  data  quantitatively  equivalent  to  those 
of  a  distribution  of  formal  potentials,  (iii)  A  kinetic  dispersion  model  based  on  a  Gaussian 
distribution  of  mnneling  distances  (or  equivalently  electronic  coupling  parameter)  from  the 
electrode  surface  is  also  evaluated.  This  model  predicts  curved  potential  step  log[i]  v^.  time 
plots,  and  in  analysis  of  log[kApp,,]  V5.  rj  plots,  undistorted  results  for  X  but  alteration  of  the 
apparent  k”. 
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INTRODUCTION 

The  extraordinary  order  of  electroactive  self-assembled  alkanethiol  monolayers 
(SAMs)  chemisorbed  on  gold  has  provoked  intense  interest^'^^  in  using  them  to  study  effects 
of  distance,  order,  and  chemical  environment  on  heterogeneous  electron-transfer  events. 
Among  the  kinetic  investigations  of  these  systems,  those  of  Chidsey'  and  Finklea^  stand  out. 
Both  groups  interpreted  log[kApp]  vs.  over-potential  (rj)  plots  with  the  Marcus^*  free  energy- 
rate  relation,  obtained  rate  constants  (k*’)  and  reorganizational  energy  barriers  (X)  to  electron 
transfer,  and  found  similar  electronic  coupling  coefficients  ((3  =  1.02  to  1.1  per  CH2)  and 
symmetrical  Tafel  plots  supporting  a  through-bond  mnneling  mechanism.  Nearly  ideal, 
kinetically-homogeneous  behavior  was  observed'-^^  in  aqueous  acid. 

Electroactive  SAMs  on  gold  in  aqueous  electrolyte  have  been  studied  over  a  limited 
temperature  range,  and  investigations  in  non-aqueous  solvents  have  encountered  monolayer 
instab ility.^'^'^®’^^  This  is  unfortunate  since  temperature  variation  is  another  route  to 
reorganization  barrier  determination,  and  being  non-adiabatic,  these  electron-transfer 
reactions  should  depend  on  dielectric  energetics  not  dynamics.  Nonaqueous  solvents  are 
more  amenable  to  measurements  at  lowered  temperatures  and  also  may  alter  monolayer  self¬ 
organization  in  insightful  ways.  Finklea  has  observed  non-ideal  cyclic  voltammetry  and  non¬ 
linear  potential  step  ln[i]  vs.  time  plots  for  [Ru(NH3)5Py]^'^''^'^  monolayers^*^,  which  gave 
relatively  stable  responses  in  several  organic  solvents.  Reorganization  energy  barriers 
determined  from  potential  step-derived  log[kApp]  vj'.  -q  plots  were  slightly  lower  in  less  polar 
solvents,  did  not  correlate  well  with  dielectric  continuum  theory,  and  were  suggested  to  be 
influenced  by  monolayer  disorder.  These  results  are  consistent  with  a  dispersion  of  kinetic 
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rates  in  the  organic  solvents. 

As  part  of  a  project  on  electron-transfer  kinetics  at  superconductor  electrodes^®,  we 
have  measured'*'^'  the  electron  transfer  kinetics  of  ferrocene-terminated  SAM’s  over  a  115K 
to  175K  temperature  range  in  a  binary  organic  cryoelectrochemical  solvent.  We,  as 
Finklea,^'^  also  observe  non-ideal  cyclic  voltammetric  and  potential-step  behavior,  and  have 
sought  to  examine  the  consequences  of  this  kinetic  dispersity  in  the  determination  of  kinetic 
parameters.  (Specifically,  the  ferrocene-terminated  monolayers  exhibit  broadened  cyclic 
voltammetric  waveshapes  and  curved  log[i]  V5.  time  plots  in  potential  step  experiments.)  We 
find  that  the  temperature  dependance  of  k®  derived  from  measured  cyclic  voltammetric 
^EpEAK  values  yields  reorganization  energy  barriers  (X)  and  (extrapolated)  k®273  rate  constants 
that  are,  for  several  SAM  alkyl  chain  lengths,  reasonably  consistent  with  dielectric  continuum 
theory  and  the  aqueous  phase  results  of  Chidsey’*^.  Analysis  of  ferrocene  SAM  voltammetric 
waveshapes  and  of  potential  step  log[kApp]  vs.  over-potential  plots,  on  the  other  hand,  gave 
much  lower  results  for  X,  and  anomalous  variations  with  temperature  and  potential  sweep 
rate. 

Hypothesizing  that  the  latter  methods  fail,  at  least  in  part,  because  of  aspects  of 
kinetic  dispersion,  we  model  effects  of  kinetic  dispersion  on  determinations  of  rate  constants, 
reorganizational  barrier  energies,  and  a  distance  transmission  coefficient.  This  paper 
presents  calculations  based  on  a  Gaussian  distribution  of  formal  potentials  of  the  redox  sites 
in  the  electroactive  SAM.  The  distribution  of  formal  potentials  is  selected  to  match  the 
broadened  voltammetric  waveshape.  The  thermodynamically-based  distribution  is  shown  to 
lead  to  an  apparent  dispersion  in  electron  transfer  rate  constants,  and  upon  analysis  by 
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log[kApp]  v^.  over-potential  plots,  low  values  of  determined  X.  We  also  consider  a  real 
kinetic  dispersion  by  modeling  a  Gaussian  distribution  of  electron-transfer  tunneling 
distances,  again  analyzing  its  effects  on  log[kApp]  vs.  over-potential  plots. 

EXPERIMENTAL 

Chemicals.  CpFeCpC02(CH2)j-SH  (Cp  =  cyclopentadienyl;  x  =  8,  12,  16)  were  prepared 
by  a  literature  method. Octanethiol,  dodecanethiol,  hexadecanethiol  (Aldrich),  HCl 
(Mallinckrodt),  HNO3  (EM  Science),  butyronitrile  (PrCN,  Aldrich,  99+%),  absolute  ethanol 
(AAPER  Alcohol  and  Chemical  Co.),  and  tetra-n-butylammonium  hexafluorophosphate 
(BU4NPF6,  Fluka,  Puriss  grade)  were  used  as  received.  Ethyl  chloride  (EtCl,  Linde)  was 
condensed  into  a  Schlenk  tube  and  stored  at  room  temperamre.  Water  was  purified  with  a 
Bamstead  NANOpure  system. 

Electrode  Fabrication.  Working  and  pseudo-reference  electrodes  were  prepared  by  soldering 
Teflon  shrouded  silver  coated  copper  wires  to  the  ends  of  a  2  mm  dia.  Ag  rod  (Johnson 
Matthey  Electronics,  99.999%)  and  0.5  mm  dia.  Au  wire,  respectively,  and  encapsulating 
these  side-by-side  in  a  cylinder  of  insulating  epoxy  (Shell  Epon  828,  m-phenylenediamine 
curing  agent,  cured  overnight  at  70°C).  The  cylinder  was  sanded  (Buehler  600  grit)  to 
expose  the  Ag  and  Au  disk  surfaces,  which  before  each  use  were  polished  with  successively 
finer  alumina  powder  (Buehler)  aqueous  slurries  down  to  0.3  /xm  grit  and  sonicated  for  five 
minutes  in  NANOpure  water. 

Chemisorbed  Monolayer  Preparation.  The  Au  electrodes  were  etched  in  dilute  aqua-regia 
(3:1:6  H2O,  HCl,  HNO3)  for  5  minutes,  rinsed  with  copious  water  then  ethanol,  immersed  in 
the  coating  solution  for  a  minimum  of  24  hours,  and  then  rinsed  thoroughly  with  ethanol. 
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Ethanolic  coating  solutions  contained  of  1:3  mole  ratios  (1  mM  total  concentration)  of 
(ferrocenylcarboxy)-octanethiol  and  n-octanethiol,  or  (ferrocenylcarboxy)-dodecanethiol  and 
dodecanethiol,  or  (ferrocenylcarboxy)-hexadecanethiol  and  hexadecanethiol.  Following 
Chidsey^’’,  co-chemisorption  with  n-alkane  thiols  has  the  purpose  of  reducing  ferrocene- 
ferrocene  interactions  and  film  defect-based  lateral  ferrocene-ferricenium  electron  transfers. 
The  Ag  pseudo-reference  was  scraped  to  expose  a  fresh  metal  surface.  No  incubation  of 
chemisorbed  films  in  diluent  solution'*’  was  employed. 

Electrochemical  Measurements.  Because  EtCl  (bp  285K)  is  a  gas  at  room  temperature  and  to 
avoid  possible  thiol  desorption,  subsequent  manipulations  were  conducted  at  dry  ice 
temperature.  The  working-reference  assembly  (vide  infra)  was  rinsed  with  chilled 
butyronitrile  and  fitted  into  a  slotted  stainless  steel  cylinder,  the  closed  Pt-faced  end  of  which 
served  as  the  counter  electrode.  The  working-counter  electrode  separation  was  adjusted  to 
ca.  0.5  mm  to  minimize  cell  resistance.  The  cell  was  filled  with  chilled  2:1  (v:v) 

EtChPrCN  electrolyte  solution  (0.075  M  in  BU4NPF6)  and  the  three  electrode  assembly 
inserted  into  the  electrochemical  cell  which  was  then  placed  into  a  liquid  helium  cryostat 
(Janis  Model  6CND-NVT,  Lakeshore  Model  805  temperature  controller)  and  brought  to  the 
desired  temperature  (in  the  present  study,  120K  to  150K,  +0.5K). 

Cyclic  voltammetry  (CV)  and  potential  step  experiments  were  carried  out  using  a 
potentiostat  of  local  construction  and  a  microcomputer  (Micro  Systems  Engineering  486)  with 
a  home  made  program  for  potential  control  and  data  acquisition.  A  home  built  active  filter^^ 
was  employed  to  improve  S/N  and  a  positive  feedback  control  circuit ^  was  employed  to 
effect  instrumental  positive  feedback  uncompensated  resistance  correction.  The 
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uncompensated  resistance  correction  could  be  estimated  from  the  setting  on  the  feedback 
compensation  module.  Actual  uncompensated  resistances  were  measured  by  AC  impedance 
in  each  experiment  at  the  temperature  of  interest  with  a  Schlumberger-Solartron  Instruments 
Model  1255  frequency  response  analyzer  and  a  Model  1286  potentiostat.  The  quantity  of 
electroactive  chemisorbed  ferrocene,  Fpc,  was  measured  from  the  charge  under  the 
ferricenium  reduction  wave  in  a  slow  potential  scan  {e.g.,  5-10  mV/s)  cyclic  voltammogram. 

Potential  steps  were  initiated  from  the  cyclic  voltammetrically-determined  surface- 
confined  ferrocene®'^  formal  potential,  E®’  and  are  thus  equal  to  the  applied  over-potential,  -q. 
The  measured  transients  were  typically  0.2-20  secs  long  depending  upon  the  temperature  and 
the  size  of  the  potential  step.  The  applied  potential  was  remrned  to  E°’  and  the  film 
equilibrated  for  2-5  min.  to  re-establish  the  original  ferrocene/ferricenium  population  before 
initiation  of  the  next  over-potential  step. 

Background  current-time  transients  were  obtained  in  each  experiment  by  carrying  out 
potential  steps  of  the  same  magnitude  and  sign  as  77,  in  the  flat,  featureless  double  layer 
potential  regions  adjacent  to  the  ferrocene  potential  region  (at  potentials  well  negative  of 
E°’)-  Currents  shown  are  thusly  background-corrected. 

Simulations.  All  simulated  voltammograms  and  current-time  transients  were  calculated  on  a 
486DX  PC  with  a  standard  spreadsheet  program  and/or  a  home  made  pascal  program'‘^ 
Simulated  voltammograms  and  potential  step  current-time  transients  for  a  Gaussian 
distribution  of  E°’  were  generated  with  21  separate  values  (enough  so  that  the  results  are 
unchanged  for  larger  numbers  of  values;  at  smaller  numbers  of  values  this  is  not  so)  of 
voltammetric  peak  potentials  and  Efwhm»  ^nd  of  rate  constants  as  a  function  of  charge 


segment  in  plots  of  log(i)  vs.  time. 

RESULTS  AND  DISCUSSION 

Cyclic  Voltammetric  Waveshape  Broadening.  A  typical  cyclic  voltammogram  at  135K  from 
our  investigation^^^  of  (ferrocenylcarboxy)dodecanethiol+dodecanethiol  mixed  monolayers 
chemisorbed  on  polycrystalline  gold  is  shown  in  Figure  1 .  The  oxidation  and  reduction 
waves  are  well-formed  and  similarly  shaped,  with  an  average  peak  width  at  half-maximum 
(Efwhm)  of  117  mV.  The  large  oxidation-reduction  peak  separation  (AEpe^k  =  310  mV)  is 
characteristic  of  slow  electron  transfer  kinetics  of  an  adsorbed  electroactive  species.  The 
ferrocene  coverage  Fp,  =  1.3x10'^'’  mol/cm^  is  ca.  1/3  of  a  full  monolayer.^^  The 
monolayers  are  exceptionally  stable  at  low  temperature,  with  Tp;.  and  voltammetry 
quantitatively  reproducible  for  as  long  as  a  week  or  more.  The  monolayers  are,  however, 
voltammetrically  unstable  at  room  temperature,  and  desorb  quickly.^ 

The  Figure  1  voltammetry  differs  from  room  temperature  aqueous  observations'®  *’  in 
that  the  double  layer  charging  current  is  larger  and  the  ferrocene  waveshape  is  broadened 
(Efwhm)-  aqueous  electrolyte,^'*  the  double-layer  capacitance  of  n-alkanethiol  monolayers 
on  Au(lll)  follows  simple  Helmholtz  behavior  and  for  a 

(ferrocenylcarboxy)dodecanethiol-l-dodecanethiol  monolayer  like  that  in  Figure  1  is^*  1.7 
AtF/cm^  The  double  layer  capacitance  for  the  monolayer  in  Figure  1,  estimated  from 
voltammetric  charging  currents  at  potentials  more  negative  than  the  ferrocene  wave,  is  18 
/xF/cm^  in  the  EtCl/PrCN  solvent,  over  10-fold  larger  (Table  I).  Such  increases  in  double¬ 
layer  capacitances  in  organic  solvents  have  been  seen  before. We  infer  that  the 
monolayer  in  Figure  1  has  significant  permeability  to  ions  and/or  solvent,  and  probably  an 
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associated  level  of  disorder.  The  capacitances  of  (ferrocenylcarboxy)- 
octanethiol + octanethiol ,  (ferrocenylcarboxy)-dodecanethiol + dodecanethiol ,  and 
(ferrocenylcarboxy)-hexadecanethiol+hexadecanethiol  monolayers  in  EtCbPrCN  all  increase 
with  increasing  temperature  (Table  I).  In  fact,  the  capacitances  of  the  C8  and  C12 
monolayers  are  both  similar  and  not  much  smaller  than  that  for  a  bare  surface,^®  the  double 
layer  capacitance  of  which  also  increase  with  increasing  temperature.^^  The  capacitance  for 
these  monolayers  seems  mainly  dominated  by  ionic  permeability  rather  than  the  absolute 
thickness  of  the  layer  anticipated  based  on  the  known  strucmre  of  similar  layers. 


i  =  nFAkJ^]^  exp[( 


RT 

-{\-a)nF 


k. 

)(-)] 


(la) 


=  kS  exp( 


(l-a)rjF. 

RT 


(lb) 


The  more  readily  analyzed  aspect  of  the  Figure  1  voltammogram  is  its  unusually  large 
EpwHM-  Calculations  show  that  a  waveshape  broadening  of  117  mV  at  135K  is  much  greater 
than  expected  from  the  degree  of  slow  electron  transfer  kinetics  indicated  by  the  over¬ 
potential  (AEpeak/2).  Calculations  using  Butler- Volmer  kinetics  for  oxidative  surface 
voltammetry  according  to  Laviron^^,  predict  that  Efwhm  should  be  only  57  mV  at  135K  and  a 
=  0.5,  for  any  value  of  v,  kb,  tj,  and  Fr*.  Other  comparisons  of  experimental  and  calculated 
Efwhm  are  shown  in  Table  I.  The  Efwhm  predicted  from  Equations  la  and  b  shows  only  a 
weak  dependence  on  T  (Table  I,  column  headed  Efwhm >  BV)  or  a.  The  experimental  Efwhm 
values  (Table  I)  increase  at  lower  temperatures  and  with  increasing  chain  length.  We  will 


discuss  these  results  in  more  detail  later  on. 
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Efwhm  for  a  surface  voltammetric  wave  can  also  be  predicted  with  the  Marcus-DOS 
kinetic  model  described  by  Tender  et  al‘‘*  and  Creager  and  Weber'*'’  for  cyclic  voltammetry, 
based  on  that  by  Chidsey  for  potential  steps*^  This  model  avoids  the  assumptions  of  Butler- 
Volmer  theory  about  the  relative  values  of  reaction  free  energy  (ij)  and  reorganization  energy 
barrier  (X),  and  accounts  for  the  continuum  of  electronic  states  in  the  metal  electrode  by 
integration  of  the  Marcus  relation^*^  '’  over  a  range  of  energies  about  the  Fermi  level  of  the 
metal  electrode, 


^oxjred,!} 


exp 


-(x±^y 


4X 


1  +  exp(jc) 


dx 


(2) 


The  outer  sphere  X  estimated  from  the  dielectric  continuum  equation^® 


e^N  1 

_ 

1  ] 

1  1 

_ 

SttCq  ' 

2d, 

\^ap 

gives  X  =  0.75  eV  at  120K  and  0.77  eV  at  150K®^  in  the  EtCl/PrCN  mixed  solvent.  (These 
X  values  are  near  the  result,  0.89  eV,  of  an  activation  analysis  of  the  (ferrocenylcarboxy)- 
dodecanethiol  monolayer  rate  constant  reported  elsewhere®'^)  Using  X  =  0.80  eV,  and 
selecting  k°  =  1.03x10®  s  *  and  /xp  =  8.54x10®  s  'eV  *  so  that  a  10  mV/s  voltammogram 
calculated  from  Equation  2  has  the  same  AEp^^K  (310  mV)  as  observed  in  Figure  1,  gives  a 
predicted  value  of  Efwhm  =  70  mV  at  135K.  The  calculated  Epwhm.M-DOS  is,  while  larger 
than  the  Efwhm>BV,  again  more  narrow  than  that  of  the  Figure  1  experiment  (as  are  other 
comparisons  between  Equation  2  and  experiment,  shown  in  Table  I).  Thus,  the  Efwhm 
predicted  from  both  kinetic  models  (Equations  1,2)  are  significantly  smaller  than  the 
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experimental  value  as  long  as  a  single  value  of  rate  constant,  X,  or  is  assumed  as  above. 
Kinetic  Dispersion  Modelled  in  Cyclic  Voltammetry  as  a  Distribution  of  Formal  Potentials. 
The  voltammetric  waveshape  broadening  in  Figure  1  and  non-linear  potential  step  ln[i]  vs. 
time  plots  {vide  infra,  Figure  4)^*  appear  to  reflect  a  dispersion  of  the  rates  at  which  the 
ferrocene  sites  react.  There  can  be  several  sources  of  such  kinetic  dispersion.  We  consider 
first  a  model  in  which  the  kinetic  dispersion  arises  from  a  spread  in  the  formal  potentials, 

E°’,  of  the  surface  ferrocene  sites.  Such  a  spread  of  thermodynamic  potentials  could  arise 
from  uneven  ion-pairing  (Fc^,  PF^)  or  ferrocene  site-site  interactions,  or  from  a  distribution 
of  monolayer  strucmres  that  induce  uneven  interactions  or  incur  a  range  of  solvation  shells. 
Whatever  its  source,  a  spread  in  E”’  will  lead  to  an  apparent  spread  or  dispersion  in  electron 
transfer  rates  and  analyzed  rate  constants,  as  we  show  in  this  section. 

Definitive  proof  of  a  distribution  of  formal  potentials,  to  distinguish  it  from  other 
potential  sources  of  kinetic  heterogeneity,  could  be  gleaned  from  a  slow  scan  rate 
voltammetric  experiment  so  as  to  give  a  reversible  response.  In  the  case  of  an  E'” 
distribution,  the  reversible  voltammetric  peak  shape  should  remain  broad,  but  it  is  difficult  at 
low  temperature  to  scan  slowly  enough  to  attain  this  condition.  However,  the  same  films 
smdied  here,  when  examined'*®  at  room  temperature  in  aqueous  electrolyte,  are  both  close  to 
reversibility  and  still  broad.  This  observation  tends  to  eliminate  other  possible  sources  of 
waveshape  broadening  such  as  double  layer  effects  as  discussed  later.  That  is,  there  is  a 
strong  possibility  that  an  E°’  distribution  exists  for  the  monolayer  results  discussed  here. 

A  kinetic  heterogeneity  based  on  a  thermodynamic  distribution  (potentially  Gaussian) 
of  energy  gaps  between  ground  and  excited  state  levels  for  donor-acceptor  molecules 
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undergoing  photochemically-induced  electron  transfer  has  been  discussed  by  Gudowska- 
Nowald®  to  explain  experimentally  observed  inhomogeneous  line  broadening  and  non¬ 
exponential  decay.  This  is  analogous  to  an  electroactive  monolayer  with  a  distribution  of  E*” 
and  thus  non- ideal  voltammetry  and  potential  step  behavior. 

We  model  a  distribution  of  E°’  with  a  standard  Gaussian  expression 

f(y)  =  - / -  exp-[(y-p)2/2(a(E:°'))2 

a(E  ^)f(2Tr)  (4) 

where  (t(E°’)  is  the  standard  deviation  of  the  distribution.  Voltammetric  curves  were 

calculated  with  the  Marcus-DOS  Equation  2  to  match  experimental  voltammograms  by 

selecting  appropriate  values  of  a(E®’)  and  summing  a  minimum  of  21  separate 

voltammograms  with  E^’’  values  lying  between  ±3cr(E‘”)  and  scaled  by  the  Gaussian  function. 

A  clarification  regarding  rate  constants  determined  with  Equation  2  from  cyclic 

voltammograms  exhibiting  a  Gaussian  distribution  of  E*”  is  useful  at  this  point.  We  consider 

two  kinds  of  rate  constants,  the  first  being  the  mid-point  average  rate  constant  k^p^,  which 

is  that  with  which  a  voltammogram  can  be  simulated  using  Equation  2  and  a  Gaussian 

distribution  of  E°’  to  fit  both  the  experimental  voltammogram ’s  shape,  Ej^hm  and  AEpeak> 

and  the  second  is  the  average  rate  constant  kV,  which  is  the  rate  constant  determined  from 

the  best  fit  of  AEpeak  values  (experimental  or  simulated)  to  values  predicted  by  Equation  2, 

assuming  a  homogeneous  population  of  reacting  sites.  The  latter  procedure  has  been 

employed'^’^^  in  our  previous  assessments  of  experimental  monolayer  kinetics,  and  involves  no 

modelling  of  the  kinetic  dispersion. 

Figure  2  shows  a  calculated  voltammogram  based  on  the  distribution  ct(E°’)  =  55 
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mV,  impmpa  =  1-48x10^  eV'^s  \  X=  0.80  eV,  T  =  135K,  and  the  mid-point  average  rate 
constant  k^MPA  =  1-78x10"'  s  '.  The  excellent  fit  to  the  experimental  voltammogram  of 
Figure  1  demonstrates  that  an  E*”  distribution  is  a  viable  explanation  for  the  observed 
voltammetric  broadening.  The  minor  differences  between  the  two  voltammograms  in  Figure 
2  appear  to  result  from  the  experimental  background  capacitance  current.  While  it  is 
possible''’^  to  fit  the  experimental  waveshapes  with  Equation  2  without  assuming  an  E°’ 
distribution,  this  requires  the  use  of  values  of  X  much  smaller  (0.40eV)''  than  that  predicted 
by  Equation  3  or  obtained  from  an  activation  plot.^‘ 

The  value  of  kV  =  1.03xl0‘^  s  '  for  a  AEp^ak  (only)  best  fit  to  Figure  1  is  5.8-fold 
larger  than  the  k^p^  determined  above  (1.78x10"*  s  ').  Thus,  in  this  example,  the  inclusion 
of  a  Gaussian  E°’  distribution  into  the  kinetic  analysis  of  experimental  data  produces  a  "mid¬ 
point  average"  rate  constant  (k^p^)  that  is  5.8-fold  smaller  than  an  "average"  analysis^' 
based  on  AEpeak  values  alone  (k°Av  )  and  assuming  no  E°’  dispersion  at  all.  We  assume  that 
the  E°’  of  the  "mid-point  average"  ferrocene  site  in  the  Gaussian  distribution  is  equal  to  the 
true  E"’  (i.e.,  for  no  dispersion,  and  where  rj  =  0  at  E  =  O.OV),  the  value  of  k^^PA  thus 
represents  the  slowest  rate  constant  of  the  kinetic  distribution  since  all  other  sites  experience 
1/  >  0  at  E  =  0  and  thus  have  larger  rates  (see  Table  II).  In  this  sense  it  is  unsurprising 
that  the  AEpE^K-only  evaluation  (k^^v  )  seems  to  moderately  favor  the  "faster"  reacting  sites 
in  a  distribution  modelled  as  a  Gaussian  spread  of  E°’  values. 

Although  the  above  difference  between  k^^PA  and  k^^v  is  not  very  large,  it  is 
important  to  examine  whether  it  is  potential  scan  rate,  temperature  or  chain  length  dependent 
since  that  would  cause  the  kinetic  dispersion  to  bias  the  measurement  of  X  (from  an 
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activation  plot)  and  then  /3  (from  a  plot  of  ln(rate  constant)  vs.  chain  length). 

That  is,  it  is  important  to  know  if  k^PA  and  kV  differ  by  a  constant  or  a  changing  factor. 
Voltammograms  simulated  based  on  the  kVp^  value  but  at  much  faster  scan  rates,  and  then 
evaluated  using  only  their  AEpeak  values,  yield  results  for  kV  (Table  II)  that  differ  from 
k°MPA  by  the  same  factor  as  that  obtained  at  0.01  V/s.  To  address  the  temperature 
dependence  of  k'^^vj  comparisons  like  that  made  at  135K  were  simulated  at  125,  150,  and 
170K;  the  results  are  shown  in  Table  II  (see  footnote  b).  As  at  T  =  135  K,  the  kV  naethod 
favors  the  faster  reacting  (at  the  average  E°’)  sites  of  the  distribution  at  these  other 
temperatures,  but  by  nearly  the  same  factor.  The  magnitude  of  the  nearly  temperature 
independent  bias  would  cause,  in  this  example,  only  -4%  error  in  X  and  +3-fold  increase  in 
fip  values  derived  for  the  C12  monolayer  from  an  activation  plot.  The  former  difference, 
and  probably  the  latter,  are  well  within  experimental  uncertainty.  Thus,  we  conclude  that  an 
E°’ -based  kinetic  dispersion  has  little  or  no  consequence  on  measurement  of  X  from  AEpeak 
values,  and  perhaps  a  modest  increase  in  the  derived  pp.  The  ca.  6-fold  difference  between 
k^^Av  and  k°MPA  can  be  expected  to  vary  with  the  a(E°’)  of  the  particular  monolayer  system. 

Comparisons  of  the  temperature  dependencies  of  experimental  Efwhm  (Table  I)  with 
simulations  involving  E"’  distribution  (Table  II)  show  that  both  Efwhm  values  decrease 
gradually  with  increasing  temperature.  The  slight  differences  between  them  are  well  within 
experiment-to-experiment  variability,  indicating  that  uIE”’)  is  nearly  independent  of 
temperature.  The  trend  in  both  is  opposite  to  the  thermal  broadening  anticipated  from 
Efwhm>BV  and  Efwhm>M-DOS  calculations  (Table  I).  Efwhm  changes  less  with  temperature 
for  the  Marcus-DOS  than  the  Butler-Volmer  predictions  (see  Table  I),  owing  to  increased 
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peak  overpotential  with  decreased  temperature  and  subsequent  peak  broadening  that  tends  to 
offset  slightly  the  thermal  narrowing.  Kinetic  dispersity  apparently  magnifies  the  latter  effect 
in  causing  the  observed  peak  narrowing  with  increased  temperature. 

The  Marcus-DOS  calculations  of  Ep^HM  in  Table  I  predict  an  increase  in  Ef^hm  with 
increased  chain  length;  this  is  caused  by  the  decreasing  rate  constant  k”  with  increased  chain 
length  and  consequent  increased  t]  relative  to  reorganization  energy.  Similarly,  a  decreased 
chain  length  should  have  the  opposite  effect. 

To  investigate  the  potential  effects  of  chain  length  on  the  difference  between  k^^PA  and 
k°Av  (and  the  size  of  a  Gaussian  distribution  of  E®’  required  to  fit  experimental  data),  a  cyclic 
voltammogram  for  a  (ferrocenylcarboxy)-octanethiol+octanethiol  mixed  monolayer  at  135  K 
(see  Table  I)  was  fit  in  the  manner  of  Figure  2.  The  best  fit  voltammogram  produced  a  mid¬ 
point  average  k^p^  =  1.10x10'^  s‘‘  and  cr(E°’)  =  47  mV,  and  the  AEpe^k  of  this 
voltammogram  produced  an  average  k\v  =  6.0x10'^.  This  distribution  size  is  only  slightly 
smaller  than  that  for  the  C 12  monolayer  (55  mV),  indicating  a  minimal  chain  length 
dependence  of  the  E®’  distribution.  The  difference  between  k'^MPA  and  k^^v  is  5.5-fold  for  the 
C8  monolayer,  similar  to  the  5.8-fold  difference  for  the  C12  monolayer  at  the  same 
temperature.  The  small  difference  between  these  factors  would  translate  into  a  /3  that  is 
smaller  by  ca.  2%  than  that  for  homogeneous  populations  for  the  two  chain  lengths.  The 
lack  of  dependence  of  the  Id’^v  analysis  on  both  potential  scan  rate  (Table  II)  and  chain  length 
is  internally  consistent  since  both  produce  changes  in  AEpeak- 

Calculations  based  on  fitting  Equation  2  with  a  Gaussian  E“’  distribution  to 
experimental  voltammograms  are  quite  tedious.  Determinations  of  a(E°')  with  Butler- Volmer 
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kinetic  Equations  la,  lb  give  values  similar  (but  somewhat  larger,  aCE^’’)  =  71  mV  for  the 
best  fit  to  Figure  1)  to  those  determined  with  the  Marcus-DOS  Equation  2,  and  are 
computationally  simpler.  (The  difference  between  Butler-Volmer  and  Marcus-DOS  wave¬ 
broadening  occurs  even  for  homogeneous  populations  of  reacting  sites.)  Since  the  Butler- 
Volmer  based  calculations  of  a(E°’)  are  simpler,  the  experimental  a(E°’)  results  in  Table  I  for 
different  chain  lengths  and  temperamres  are  derived  in  this  manner,  using  Equations  la,  lb. 
Other  Wave-Broadening  Effects.  Other  possible  sources  of  voltammetric  waveshape 
broadening,  include  (i)  double-layer  effects,  (ii)  uncompensated  resistance,  and  (iii)  a  real 
dispersion  of  rate  constants  (vide  infra).  Double-layer  effects  for  redox-active  monolayers 
have  been  described  theoretically  by  Smith  and  White'^*  and  experimentally  observed  by 
Rowe  and  Creager^^  In  the  present  experiments  the  supporting  electrolyte  concentration  is 
relatively  small  (0.075M)  and  the  actual  ion  population  probably  even  smaller  due  to  ion 
pairing  at  low  temperatures.  Also,  Epc  is  moderately  high  (1/3  of  a  full  monolayer).  Since 
such  factors  have  in  double  layer  simulations  been  shown  to  cause  wave  broadening,  we 
calculated  voltammograms  to  explore  them."*^  The  calculations  predict  Efwhm  values  (70  mV 
at  125K  and  83  mV  at  150K)  much  smaller  than  those  observed  (Table  I),  and  additionally  a 
skewing  of  the  waveshape®®’''^  These  observations  suggest  that  double-layer  effects  are  not 
the  dominant  waveshape-broadener  in  our  experiments.  Additionally,  ion-pairing  of  surface- 
confined  ferricenium  with  PFg'  is  highly  probable  at  low  temperature;  double-layer  effects 
would  vanish  on  a  thus-neutralized  surface.  It  is  known,  for  example,  that  PFg'  ion-pairs 
with  electrogenerated  ferricenium  in  monolayers  containing  ferrocenylalkanethiols  in  aqueous 
electrolyte  at  room  temperature.^*'  Lastly,  penetration  of  ions  into  the  nominally  pure 
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hydrocarbon  region  of  the  monolayer,  as  suggested  by  the  capacitance  data  in  Table  I,  may 
screen  the  electric  charge  from  the  redox  centers  thus  eliminating  any  double-layer  effect. 

We  next  consider  uncompensated  resistance  (Runc)  a  source  of  voltammetric 
broadening.  Positive  feedback  compensation  is  utilized  in  these  experiments  to  correct  for 
uncompensated  resistance,  and  it  can  be  expected  that  the  instrumental  compensation  is 
always  less  than  100%  of  Runc  ”*^  Typical  values  of  Ry^c  that  are  corrected  for  with  the 
positive  feedback  system  are  between  1x10^  to  1x10^  ohms  in  the  temperature  range  of  125K 
to  150K.  Direct  AC  impedance  measurements  of  Rung  generally  agree  with  the  instrumental 
settings.  It  is  important  to  realize  that  even  at  slow  potential  scan  rates,  the  slow  electron 
transfer  rates  characteristic  of  low  temperatures  (and  tunnelling  barriers)  can  give  large 
AEpEAK  values  (Table  I)  and  small  peak  currents  in  voltammograms  (a  few  nanoamps  or  less), 
so  that  iRuNc  is  small  relative  to  AEpe^k  even  though  the  Rung  remaining  after  feedback 
compensation  is  by  usual  electrochemical  standards  rather  large.  As  a  result,  experimental 
voltammograms  are  often  quite  similar  with  and  without  positive  feedback.  Calculations  of 
cyclic  voltammograms  with  Marcus-DOS  kinetics  (Equation  2)  based  on  the  best  fit 
parameters  to  the  data  in  Figure  1  and  including  the  effect  of  an  Rung  2  predict 

no  additional  waveshape  broadening  (beyond  that  reported  in  Table  I)  at  a  potential  scan  rate 
of  5  mV/s,  and  an  increased  broadening  of  5  mV  at  a  potential  scan  rate  of  10  mV/s.  It  is 
therefore  unlikely  that  the  AEpeak  or  Ep^HM  results  reported  in  Table  I  are  caused  by 
uncompensated  resistance. 

Lastly,  waveshape  broadening  can  result  from  a  real  distribution  of  standard  rate 
constants,  such  as  caused  by  a  distribution  of  tunneling  distances  or  of  reorganization  barrier 
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energies;  this  is  discussed  later. 

Effects  of  a  Gaussian  Distribution  of  Ef’  on  Potential  Step  Experiments.  In  potential  step 
experiments,  a  Gaussian  distribution  of  E®’  values  provokes  an  apparent  distribution  of 
electron  transfer  rates  because  the  applied  overpotential  77  is  reckoned  with  respect  to  the 
average  E°’  of  the  redox  couple.  The  effective  overpotential  experienced  by  any  given  redox 
site  will  vary  with  its  microscopic  formal  potential,  resulting  in  a  distribution  of  reaction 
rates  over  the  ensemble  of  surface  redox  sites  and  a  consequent  decay  of  current  with  time 
which  is  not  strictly  first  order.  This  apparent  rate  dispersion  may  of  course  occur 
concurrently  with  a  real  dispersion  of  standard  rate  constants  for  different  reacting  sites. 

The  first  order  reaction  of  a  surface-confined  species  is  expected^^  to  display  an 
exponentially  decaying  current. 


i{f)  ^APP,t]Q  ^APP,if^ 


(5) 


where  k^pp,,  (s'^)  is  the  sum  of  the  forward  and  reverse  electron  transfer  rate  constants  (ko,^ , 
h  ^red.,,)  at  the  effective  overpotential  77,  and  Q  is  the  charge  passed  in  the  electron  transfer 
reaction. 


Current-time  responses  to  applications  of  effective  overpotential  (applied  potential 
relative  to  the  average  E°’)  and  for  given  values  of  and  were  calculated  following 
Equation  5.  In  these  calculations,  the  rate  constant  k^pp,,  depends  on  the  effective 
overpotential  as  in  Equation  2,  while  np  and  X  are  assumed  to  be  constant.  To  calculate  the 
current-time  curves,  the  total  charge  for  the  reacting  monolayer  was  divided  into  21 
segments,  each  with  its  own  E°’.  The  size  of  each  charge  segment  is  determined  by 
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multiplying  the  total  charge  by  the  normalized  intensity  of  the  Gaussian  function  at  each 
given  The  objective  was  to  generate  plots  of  log[kApp,J  ys.  effective  rj,  to  be  analyzed 

by  best  fits  to  Equation  2  to  see  whether  the  determined  kinetic  parameters  X,  /xp,  and  k“ 
differ  from  the  values  pertinent  to  a  surface  population  having  a  single  E°’  value  equal  to  the 
average  E°’. 

Figure  3A  and  3B  show  calculated  125K  current-time  responses,  plotted  as  log[i]  V5. 
time  according  to  Equation  5,  for  potential  steps  to  rj  =  0.0  V  and  0.25  V,  respectively,  with 
respect  to  the  average  E°’,  each  assuming  Fp^  =  1.3x10''°  mol/cm^,  and  for  a(E°’)  values 
ranging  from  0  to  100  mV.  k^pp,,  for  each  rj  was  calculated  from  Equation  2  using  X  = 

0.80  eV  and  /xp  =  1  xlO^  s'‘  eV'*.  Figure  3A(curve  E)  and  3B(curve  D)  correspond  to  no 
dispersion  in  E°’,  and  yield  linear  log[i]  vs.  time  plots.  All  others  are  curved.  The  extent  of 
curvature  increases  with  increasing  formal  potential  dispersion  (t(E°’),  i.e.,  an  increasing 
dispersion  in  the  apparent  reaction  rate.  The  difference  in  rate  (slope)  for  a(E°’)  =  60  mV 
between  the  shortest  and  longest  times  is  about  10-fold.  Clearly  the  sites  that  experience  a 
greater  effective  overpotential  generate  more  current  at  shorter  times  and  dominate  the 
response  there,  whereas  sites  experiencing  a  smaller  effective  overpotential  dominate  at 
longer  times.  The  faster  sites  dominate  over  the  time  range  shown  for  17  =  0  and  a(E°’)  > 

0,  and  the  slope  lies  above  the  actual  rate  except  at  the  largest  times,  whereas  for  17  =  0.25 
V,  the  average  slope  (ka^e)  passes  from  higher  to  lower  than  the  actual  rate  near  the  middle  of 
the  plot. 

Figure  4  shows,  for  comparison,  a  log[i]  vs.  time  plot  of  typical  experimental  data  at 
125K  and  77  =  0.2  V  for  a  monolayer  of  (ferrocenylcarboxy)-dodecanethiol  coadsorbed  with 
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dodecanethiol.  The  plot  in  Figure  4  is  curved  with  a  4-fold  difference  in  rate  from  longest  to 
shortest  times.  Other  potential  step  experiments  in  EtCl/  PrCN  at  low  temperature  have 
given  similarly  curved  log[i]  vs.  time  responses^*  regardless  of  chain  length,  temperature,  or 
applied  overpotential.  Others  have  seen  such  curvamre.^*"  The  calculations  shown  in  Figure 
3  demonstrate  that  a  distribution  of  formal  potentials  a(E°’)  selected  to  fit  a  broadened  cyclic 
voltammogram  (i.e..  Figure  2),  produces  log[i]  vs.  time  plots  that  are  curved  in  a  manner 
similar  to  that  observed  experimentally  (Figure  4). 

The  more  significant  issue  is  whether  the  curved  log[i]  vs.  time  plots  described  above 
translate  into  any  substantial  bias  when  they  are  analyzed  according  to  Equation  2  for 
apparent  values  of  k“,  X,  and  ixp.  Figure  5  shows  three  log[kApp  ,]  v^.  rj  plots  generated  from 
calculated  125K  current-time  plots  with  {t(E®’)  =  45  mV.  k^pp,,  for  each  -q  was  calculated 
from  Equation  2  using  values  of  X  and  pp  of  0.80  eV  and  1  x  10®  respectively.  The 
overpotential  range  shown  in  Figure  5  is  typical  of  that  currently  accessible  in  our 
experiments^*.  Due  to  the  curvature  of  the  log[i]  vs.  time  plots,  the  total  faradaic  reaction 
charge  was  divided  into  ten  equal  segments  and  the  rate  constant  evaluated  for  each  segment 
from  the  average  slope  during  the  time  interval  over  which  that  charge  segment  is  passed. 

This  procedure  allows  the  determination  of  rates  for  the  same  population  of  ferrocene  sites  as 
a  function  of  applied  overpotential.  Figure  5  (□,  V,  O)  correspond  to  the  indicated 
ferrocene  charge  segments,  i.e.,  the  first  0-10%  of  reacting  sites,  the  30-40%  of  reacting 
sites,  etc. 

Figure  5  also  shows  a  plot  (O)  for  an  ideal  homogeneous  population  (a(E°’)  =  0  mV) 
of  redox  sites.  This  ideal  curve  (O)  lies  below  the  calculated  ones  and  exhibits  less  curvature 
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at  higher  applied  rj  values.  The  calculated  log[kApp,,]  V5.  rj  plot  for  each  charge  segment  in 
Figure  5  was  compared  to  Equation  2  to  extract  a  best  fit  for  k“,  X,  and  /xp.  The  results  at 
125K  are  summarized  in  Table  III.  The  best  fit  to  the  first  (0-10%  reacting)  charge  segment 
produces  X  =  0.55  eV,  pp  =  5.4x10^*  s  'eV  ’,  and  k®  =  5.0x10'^  s  ^  The  apparent  values  of 
X  and  pp  are  smaller  and  that  for  k°  larger  than  values  for  a  kinetically  homogeneous 
population  (a(E^’)  =  0  mV).  The  later  charge  segments  give  values  of  X,  pip,  and  k'’  that  are 
more  similar  to  but  still  considerably  different  from  the  ideal  values,  especially  pp.  Similar 
calculations  (Table  III)  show  that  the  distortion  in  the  derived  kinetic  parameters  is  less  at 
higher  temperature  and  substantially  greater  for  larger  alE*^’)-  It  is  significant,  however,  that 
the  average  rate  constants  for  later  charge  segments,  i.e.,  in  the  second  one-half  of  the 
reacting  population,  dijfer  from  the  ideal  standard  rate  constant  kf  by  less  than  two-fold.  For 
this  reason  we  regard  k°  data  derived  from  experimental  potential  step  results  that  exhibit 
kinetic  dispersity  (i.e..  Figure  4)  in  the  manner  of  Figure  5  as  reliable  to  a  ca.  a  factor  of 
two,  as  long  as  the  analysis  is  confined  to  the  later-reacting  surface  population.  Similarly- 
derived  values  of  X  and  pp  are  more  suspect. 

A  comparison  of  the  homogeneous  rate  and  that  derived  at  each  charge  segment  for 
simulated  potential  step  experiments  is  analogous  to  comparing  k%pA  and  k\v  derived  from 
simulated  cyclic  voltammograms.  Recall  that  best  fit  voltammograms  are  such  that  the 
difference  between  k^^PA  and  k^^y  is  approximately  6-fold;  a  similar  difference  in  rates  for 
potential  step  simulations  is  found  between  the  rate  determined  from  the  third  to  fourth 
segment  and  the  average  rate  (i.e.  the  rate  determined  with  no  distribution).  The  significance 
of  this  observation  is  tenuous,  however,  because  the  rates  determined  for  charge  segments 
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are  complicated  by  the  unusually  low  value  of  X  necessary  to  fit  the  Tafel  plot. 

Finally,  another  possible  reason  for  curved  log[i]  vs  time  plots  is  charge  effects  within 
the  electrical  double  layer  as  described  by  Weber  and  Creager^^.  As  previously  stated,  we 
consider  that  double-layer  effects  are  probably  minimal  in  our  experiments  due  to  ion¬ 
pairing.  Also,  log[i]  vs.  time  plots  exhibit  the  same  magnimde  of  curvature  for  oxidative  and 
reductive  overpotential  steps  (data  not  shown),  which  is  not  predicted  by  the  double-layer 
model. 

Kinetic  Dispersion  Modelled  as  a  Gaussian  Distribution  of  Reorganizational  Energy  Barriers 
or  Tunneling  Distances.  Kinetic  dispersion  can  of  course  be  real  and  not  apparent  as  in  the 
E°’  distribution  discussed  above.  Of  the  several  forms  of  kinetic  heterogeneity  in 
electroactive  SAM’s,  we  consider  (i)  a  distribution  of  reorganizational  barrier  energies  X 
(such  as  might  result  from  variable  ferrocene/monolayer/solution  interactions),  and  (ii)  a 
distribution  of  tunnelling  distances  or  electronic  coupling  factors  (such  as  might  arise  from 
monolayer  disorder). 

Consider  a  distribution  in  X.  A  monolayer  with  ferrocene  sites  that  exhibit  a 
distribution  of  X  values  (but  a  uniform  E*^’)  will  generate  a  kinetic  dispersion  that  is  real  but 
which  in  appearance  is  very  similar  to  the  apparent  dispersion  caused  by  a  distribution  in  E'^’ 
at  constant  X.  The  similarity  occurs  because  E°’  (as  rj)  and  X  both  appear  in  the  exponential 
of  Equation  2.  E°’  and  X  dispersions  under  conditions  approaching  reversibility  {i.e.,  very 

slow  potential  scan  rates,  higher  temperatures,  shorter  alkyl  chains)  would  however  produce 
differing  voltammetry;  for  the  former  dispersion,  broadened  waveshapes  would  persist, 
whereas  for  the  latter  the  waveshape  would  become  ideal.  It  seems  likely  that  in  reality  X  or 
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E°’  dispersity  will  occur  together,  both  contributing  to  the  kinetic  behavior,  and  making  their 
deconvolution  difficult. 

Kinetic  dispersion  can  also  appear  as  a  distribution  in  the  electronic  coupling 
parameter  /i,  which  depends^*  exponentially  on  distance,  d. 


p  =  Po  Q^i-pid-do)) 


(6) 


where  po  represents  the  coupling^**’  at  the  smallest  electrode-redox  molecule  separation 
distance  d^,  and  /3  is  the  electronic  coupling  parameter  or  distance  decay  constant.  Thus  a 
distribution  of  either  the  distance  over  which  electron  transfer  occurs  or  the  coupling  term  /3 
(by  chain  conformational  changes'*’)  would  result  in  a  distribution  of  p  and  thus  a  dispersion 
in  the  reaction  rate.  This  dispersion  was  modelled  by  assuming  that  at  the  moment  of 
electron  transfer  a  Gaussian  distribution  of  distances  exists  over  which  electron  transfer  can 
occur.  (The  spread  in  this  distribution  is  implied  to  depend  on  the  degree  of  monolayer 
disorder.) 

Figure  6  shows  a  log[i]  vs.  time  plot  for  a  Gaussian  distribution  of  distances,  leading 
to  a  distribution  of  p  as  described  by  Equation  6,  for  a  monolayer  with  other  parameters  as 
in  Figure  1,  an  average  ferrocene  distance  of  14 A  from  the  electrode  surface,  and  an 
overpotential  of  0.25V  at  125K.  The  standard  deviation  of  the  distance,  a(d)  =  1.54A  is 
taken  as  the  length  of  a  carbon-carbon  spacing  in  an  n-alkane  and  a  constant  value  of  /3  = 
l.OSA  *.  The  simulated  transient  (Figure  6)  was  determined  with  Equation  4,  Tpc  =  1.3  x 
lO  ***  mol/cm’,  and  a  rate  constant  that  in  Equation  2  corresponds  to  X  =  O.SOeV,  and  pp^ve 
=  6  X  10®  s  *eV'*  at  an  average  distance  of  14A.  The  value  of  p  varies  from  ±3*ff(d). 
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according  to  Equation  6.  Figure  6  shows  a  pronounced  curvature,  indicative  of  kinetic 
dispersion  as  anticipated.  The  rate  determined  from  the  slope  of  Figure  5  is  greatest  at 
shorter  times  and  varies  by  over  an  order  of  magnimde. 

Tafel  plots  derived  from  a  charge  segment  analysis  of  Figure  6  (similar  to  that  of 
Figure  5)  are  shown  as  Figure  7,  with  best  fit  parameters  for  charge  segments  1,  4,  and  7 
listed  in  Table  IV.  As  expected,  the  results  vaiy  with  the  charge  segment.  Segment  4  yields 
values  of  X,  np,  and  k®  close  to  the  original  (no  dispersion)  values.  Interestingly,  X  varies 
only  mildly  with  charge  segment,  which  is  a  significant  difference  between  this  type  of 
kinetic  dispersion  and  that  due  to  a  Gaussian  distribution  of  E°’ .  The  minor  variation  in  X 
seems  reasonable  since  the  dispersion  is  in  the  pre-exponential  term,  which  contains  only  a 
weak  dependence^**’  on  X.  The  other  variables,  pp  and  k*’,  do  however  change  considerably 
from  one  charge  segment  to  the  next. 

It  is  instructive  to  examine  the  physical  aspects  of  a  self-assembled  monolayer  that 
could  lead  to  a  kinetic  dispersion  through  Equation  6.  Several  studies  have  deduced  similar 
values,  1.0  to  1.1 /CH2  for  electron  transfers  through  n-alkanethiol  monolayers  on 
gold.*®  **’’^®  '’*  ***  (3  was  also  determined  to  be  independent  of  the  electrode  potential,  consistent 
with  a  through-bond  tuimeling  mechanism.  Similar  values  result  from  kinetic  studies 
with  immobilized*^’*"’^  and  with  soluble  redox  probes'**’'*^  at  SAM’s,  suggesting  that  the 
mechanism  for  electron  transfer  across  the  monolayer  is  similar  for  these  different  physical 
situations,  i.e.,  the  tunnelling  barriers  are  similar  for  electron  transfers  through  coadsorbed 
n-alkanethiols  and  through  the  chain  of  the  redox  probe.  In  fact,  recent  results  by  Finklea^^^ 
support  the  notion  of  electron  transfer  through  the  diluent  alkanethiol  in  mixed 
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alkanethiol/redox-active  alkanethiol  monolayers. 

Since  the  termini  of  (ideally,  all  trans)  SAM  alkane  chains  are  known  to  have  gauche 
defects and  defects  there  and  in  the  chain  interior  may  be  more  pronounced  in  organic 
solvents,  it  is  reasonable  to  suppose  that  variances  can  occur  in  distance  d  or  in  the 
orientation  of  the  chains  (i.e. ,  0).  Super-exchange  calculations'^^  for  the  electron  transfer  of 
donor-acceptor  molecules  predict  substantial  changes  in  coupling  with  alkyl  orientation.  For 
example,  a  single  gauche  bond  in  the  middle  of  an  all-trans  n-alkane  chain  is  calculated'*^’’  to 
decrease  the  coupling  matrix  (Hab)  by  a  factor  of  twenty  (from  76-1591cm  ’),  which 
translates  into  a  factor  of  two-fold  change  in  jS.  Gauche  defects  concurrently‘’’j  produce  a 
change  in  d  (measured  normal  to  the  surface),  so  dispersion  through  Equation  6  requires 
inter-related  a(d)  and  a(/3)  dispersions.  Experimentally,  an  alkylsiloxide  SAM  has  been 
observed'*’’’  to  exhibit  a  larger  tunneling  barrier  upon  annealing  to  introduce  gauche  defects. 
Correlating  monolayer  structure  (disorder)  and  coupling  and  thus  electron  transfer  rate 
dispersion  in  SAM’s  is  an  interesting  experimental  prospect  which  will  require  application  of 
other  structure-sensitive  techniques. 

A  more  detailed  examination  of  the  consequences  of  a  distribution  of  tunneling 
distances  is  possible,  particularly  in  regards  to  fitting  experimental  cyclic  voltammograms, 
but  since  we  have  little  structural  information  regarding  our  system,  and  since  simulated 
Tafel  plots  with  such  distributions  did  not  agree  well  with  our  experimental  results,  we  did 
not  pursue  this  notion.  It  is  clear  from  the  above  discussion  how  one  would  go  about  doing 
so,  and  may  be  more  appropriate  for  other  systems. 

Summary  A  Gaussian  distribution  of  E°’  with  a(E°’)  of  approximately  55  mV  is  consistent 
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with  the  experimentally  observed  waveshape  broadening  in  cyclic  voltammetric  experiments 
and  the  curvature  in  log(i)  y^  time  plots  derived  from  potential  step  experiments. 
Interestingly,  fitting  kinetic  parameters  to  Tafel  plots  derived  from  calculated  potential  step 
experiments  that  account  for  a  Gaussian  distribution  of  E'’’  gives  low  values  of  X  and  [xp  and 
high  values  of  k°  as  compared  to  the  true  values  (i.e.  those  expected  with  a  homogeneous 
population  of  kinetic  sites)  at  the  earliest  charge  segments,  and  approach  but  do  not  equal  the 
true  values  at  later  ones.  Thus,  analyzing  Tafel  data  derived  from  potential  step  experiments 
with  monolayers  having  kinetic  dispersion  caused  by  a  Gaussian  distribution  of  E°’  can  yield 
erroneously  low  values  of  X  and  pp  and  high  values  of  k°.  Error  in  k°  can  be  reduced  to  cg^ 
two-fold  by  reliance  on  later  charge  segments.  Results  from  cyclic  voltammetric  simulations 
give  an  average  rate  constant  roughly  6-fold  larger  than  the  mid-point  average  rate  constant 
incorporating  a  Gaussian  E°’  distribution.  This  bias  is  nearly  temperature  and  chain  length 
tindependent  and  does  not  significantly  effect  values  of  X  and  pp  determined  from  activation 
plots  or  determination  of  /?.  Finally,  a  Gaussian  distribution  of  reorganizaiton  energies, 
which  is  real  Kinetic  dispersion,  has  consequences  on  the  appearance  and  the  analysis  of  data 
qualitatively  equivalent  to  those  of  a  distribution  of  formed  potentials. 

Kinetic  heterogeneity  caused  by  a  Gaussian  distribution  of  tunneling  distances  has  also 
been  considered.  Tafel  plots  derived  from  these  simulations  result  in  best  fit  values  of  pp 
and  k°  that  are  higher  than  the  true  values  at  early  charge  segments,  nearly  the  same  at 
intermediate  segments,  and  lower  at  later  ones.  The  value  of  X  is  nearly  the  same  as  the  true 
value  at  all  charge  segments,  which  is  a  significant  difference  between  this  type  of  kinetic 
heterogeneity  and  that  of  a  Gaussian  distribution  of  E°’ .  Therefore  this  type  of  kinetic 
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heterogeneity  can,  in  principle,  be  distinguished  from  a  Gaussian  distribution  of  E°’  based  on 
its  ability  to  render  values  of  X  that  are  in  accord  with  less  dispersion-sensitive 
approaches^  '*’^^  like  cyclic  voltammetric  AEpe^k  analysis  or  activation  plots. 

Predictions  made  with  a  Gaussian  E®’  distribution  appear  to  give  the  most  consistent 
agreement  with  experimental  cyclic  voltammetric  and  potential  step  results  and  is  therefore 
favored  as  the  likely  cause  of  apparent  kinetic  heterogeneity. 
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TABLE  I.  Summary  of  Cyclic  Voltammetric  Results  as  a  Function  of  Chain  length  and  Temperature 

Monolayer  T  (K)  AEp(mV)  Epwhm  EFWHM  Epwhm  M-  (mV/s)  Q, 

Expt.  (mV)”  BV(mV)‘^  DOS(mV)‘*  (/xF/cm^) 

(mV)*  Fitted  Calc.  Calc. 


Tj- 

m 
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TABLE  II.  Best  Fit  Kinetic  Results  for  Simulated  Cyclic  Voltammetric  Waves  Exhibiting  a  Gaussian 
Distribution  of  E°’,  as  a  Function  of  Temperature.^ 


T(K) 

V  (V/s) 

kV  (s-‘)^ 

EpWHM 

125*’ 

0.01 

2.80x10-" 

118 

100 

3.20x10"* 

202 

ISS" 

0.01 

1.03x10’ 

117 

0.1 

1.08x10-’ 

128 

10 

1.13x10’ 

166 

100 

1.24x10’ 

202 

150*’ 

0.01 

6.0x10’ 

115 

100 

6.6x10’ 

198 

170^’ 

0.01 

4.15x10-’ 

109 

100 

4.0x10-’ 

194 

^  Average  rates  constants  determined  from  the  best  fit  to  the  AEpe^k  with  Equation  2  and  X=0.80  eV  for  a  simulated 
voltammogram  with  a  Gaussian  distribution  of  E®’  (ff(E®’)=55  mV,  -48x10^  eV's’,  X=0.80eV).  b  Mid-point 

average  rate  of  the  distribution,  =  1.78xl0^s  '  determined  at  135K  by  comparison  to  experiment  as  in  Figure  2.  Mid¬ 

point  average  rate  constants  at  other  temperatures,  using  pPmpa  =  1.48xl0^eV  ’s''  and  X  =  0.80eV,  are  4.18xl0'’s  '  (125K), 
1.10xl0‘^s  ‘  (150K),  7.7xlO  V'  (170K).  c  EpwHM  are  for  simulations  that  include  the  Gaussian  E°’  distribution. 


TABLE  Ill/  Best  Fits  to  Calculated  Tafel  Plots  Derived  from  a  Gaussian  Distribution  of  E 


Segment  (eV) 


fjip  (eV  ’s' 
xlO"* 


Charge 

Segmen 

f 


X 

(eV) 

pp  (eV'^s" 

xlO-^ 

kVs- 

xlO' 

0.50 

5.1 

1.5 

0.55 

5.95 

5.5 

0.61 

9.53 

2.2 

0.65 

13.1 

1.2 

0.67 

13 

0.75 

0.69 

11.1 

0.4 

0.71 

11 

0.25 

0.8 

100 

0.28 

Charge 


pp  (eV'^s'  k®  ($■’) 


Segment  (eV)  ‘) 

'  X  10-^ 


0.55 


0.57  4.8 


0.63  10.1 


0.66  14.4 


0.71  28.4 


0.73  34.8 


0.75  38.3 


0.80  100 


Best  fit  to  simulations  with  Equation  2.  Total  charge  divided  into  10  equally  sized  segments.  All  current-time  transients 
simulated  to  generate  the  values  in  the  table  were  determined  with  Equation  2,  4,  and  5  with  X  =  0.80  eV  and  fip  =  1.0 
10^  eV  's-'. 

T  =  125K,  ff(E“’)  =  45  mV 
T  =  125  K,  aCE”’)  =  60  mV 
'T  =  150  K,  ffCEO’)  =  35  mV 
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TABLE  IV.“  Best  Fit  Kinetic  Parameters  to  Tafel  Plots  Derived  from  a  Gaussian  Distribution  of  Redox 
Probe  Distances 


Segment 

fip  (eV-'s-‘) 
xlO" 

X(eV) 

1 

9.4 

0.75 

8.5 

2 

4.4 

0.75 

4 

3 

5.1 

0.78 

2.3 

4 

3.5 

0.78 

1.6 

5 

2.5 

0.78 

1.1 

6 

1.6 

0.78 

0.7 

7 

1 

0.78 

0.45 

8 

0.56 

0.78 

0.25 

No 

Distribution 

6 

0.80 

1.7 

^  Best  fit  to  simulations  with  Equation  2.  Total  charge  divided  into  10  equally  sized  segments.  All 
current-time  transients  simulated  to  generate  the  values  in  the  table  were  determined  with  Equations 
2,4,5  and  7  with  a(d)=  1.54A,  X  =  0.80  eV,  fip^,  =  6.0  x  10"  eV  ’s  ',  T  =  125K,/3  =  l.OSA  ',  and 
d^^e  =  14A. 
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FIGURE  LEGENDS 

Figure  1.  Cyclic  voltammetry  of  a  CpFeCpC02(CH2)i2SH/CH3(CH2)iiSH  mixed  monolayer  on 
2.0x10"^  cm^  polycrystalline  gold  in  2:1  EtCl:PrCN  with  0.075M  BU4NPF6  (prepared  from  a  3:1 
CH3(CH2)iiSH/CpFeCpC02(CH2)i2SH  ethanolic  solution  at  1  mM  total  thiol  concentration).  Fpc 
=  1.3x10'^°  mol/cm^. 

Figure  2.  Cyclic  voltammogram  (anodic  wave  only)  from  Figure  1  (open  circles),  and  a 
calculated  voltammogram  (open  squares,  using  Equations  2  and  4  with  an  E”’  distribution  of 
a(E“’)  =  55  mV,  k^pA  =  1.78x10-^,  mpmpa  =  1.48x10^  eV'^s  ',  X  =  0.80  eV,  j'  =  10  mV/s,  T 
=  135K,  Fr*  =  1.3xl0‘^'’  mol/cm^;  21  separate  voltammograms  were  scaled  by  the  Gaussian 
function  and  added  together  to  generate  a  smooth  simulated  voltammogram.) 

Currents  are  normalized  to  ip. 

Figure  3.  Log[i]  time  plots  generated  from  simulated  potential  step  experiments  with  several 
distributions  of  E*”  values  (e.g.  (7(E°’))  following  Equations  4  and  5.  The  potential  step  r]  is 
equal  to  the  average  E°’  minus  the  applied  potential.  The  rates  for  each  transient  were  calculated 
from  Equation  2  with  X  =  0.80  eV  and  /xp  =  1.0  x  10^  eV  's  '.  The  quantity  of  ferrocene 
consumed  in  the  reaction  isT^,  =  1.3  x  lO"*'’  mol/cml  (A)  tj  =  O.OV  at  T  =  125K,  Curves  A- 
E  are  ff(E°’)  =  100,  80,  60,  40,  0  mV.  (B)  rj  =  0.25V  at  T  =  125K,  Curves  A-D  are 
ff(E°’)  =  100,  60,  40,  OmV. 

Figure  4.  (A)  Experimental  current  time  transient  at  T=125K  and  rj=0.20  V  for  a 
CpFeCpC02(CH2)i2SH/CH3(CH2)iiSH  mixed  monolayer  on  polycrystalline  gold  in  2:1 
EtChPrCN  with  0.075M  BU4NPF6  (prepared  from  a  3:1  CH3(CH2)„-SH/CpFeCpC02(CH2)i2SH 
ethanolic  solution  at  1  mM  total  thiol  concentration).  7.95x10^  Q  positive  feedback 
compensation  and  background  current-subtracted.  (B)  Log(i)  time  plot  of  the  data  in  (A). 
Figure  5.  Log(kApp,„)  y^  plots  generated  from  calculated  current  time  transients  with  an  E®’ 
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distribution,  aCE”’)  =  45  mV,  following  Equations  2,4,  and  5.  k^pp,,  calculated  at  each 
overpotential  with  T  =  125K,  X  =  O.SOeV  and  /xp  =  IxlO^eV'^s  ^  The  quantity  of  ferrocene 
consumed  in  the  reaction  is  Fpc  =  1.3x10'^*^  mol/cm^.  Charge  segments  1  (squares),  4 
(triangles),  and  7  (diamonds)  correspond  to  0-10,  30-40,  60-70%  increments  of  reacting 
population.  Best  fits  to  each  segment  are  shown  as  lines  through  the  calculated  rates:  segment 
1,  X  =  0.55eV,  fxp  =  5.41xl0'‘  eV'^s  ’,  and  k°  =  5.0x10'^  s'^;  segment  4,  X  =  0.67eV,  pp  = 
1.74x10^  eV'^s  *,  and  k°  =  1.0x10'^  s'\  and  segment  7,  X  =  0.72eV,  pp  =  1.88xl(F  eV'^s  \ 
and  k”  =  3.4x10'^  s  ^  Ideal  curve,  a(E°’)  =  0,  shown  as  (O). 

Figure  6.  Log[i]  time  plots  from  calculated  potential  step  experiments  (rj  =  0.25V  at  T  = 
125K)  with  a  1.54A  distribution  of  redox  probe  distances  from  the  electrode  surface,  following 
Equations  2,  4  and  7.  The  rates  for  each  transient  were  determined  from  Equation  2  with  X  = 
0.80  eV,  ppjve  =  6.0x10®  eV''s  ‘,  d^ve  =  14A,  <T(d)  =  1.54A,  /3  =  1.08A''.  The  quantity  of 
ferrocene  consumed  in  the  reaction  is  Tpc  =  1.3x10  '°  mol/cm^. 

Figure  7.  Log(kApp,,)  v  plots  from  current  time  transients  calculated  for  a  1.54A  distribution 
of  d,  following  Equations  4  and  7,  for  three  charge  segments,  k^pp,,  calculated  at  each 
overpotential  with  T  =  125  K,  X  =  0.80eV,  pp^^^  =  6.0x10®  eV  's  ',  (7(d)=  1.54A,  d^ve  =  14A, 
and  /3  =  1.08 A  '.  The  quantity  of  ferrocene  consumed  in  the  reaction  is  Tp^  =  1.3  x  10  '° 
mol/cm^.  Best  fits  to  each  segment  are  shown  as  lines  through  the  calculated  rates;  for  segment 
1,  X  =  0.75eV,  pp  =  9.43x10®  eV  's  ',  and  k°  =  8.5x10'^  s  ';  segment  4,  X  =  0.78eV,  pp  = 
3.5x10®  eV  's  ',  and  k°  =  1.6x10'^  s  ';  and  segment  7,  X  =  0.78eV,  pp  =  1.0x10® eV  's  ',  and  k° 
=  4.5x10'^  s  '.  Ideal  curve,  a(d)  =  0,  shown  as  (O). 
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